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ABSTRACT 


An  effect  of  the  atmospheric  channel  on  near-ground  propagation  of  acoustic  waves  is  considered.  The  structure 
of  the  amplitude  and  phase  fluctuations  of  acoustic  waves  propagating  through  the  turbulent  atmosphere  near  the 
ground  has  been  analyzed.  The  structure  of  the  temperature  and  wind  velocity  fields  and  vertical  profiles  of  the 
outer  scale  of  atmospheric  turbulence  have  been  experimental  lv  investigated. 
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1.  INTRODUCTION 


Investigations  into  the  structure  of  the  atmospheric  boundary  layer  and  acoustic  waves  propagating  in  it  are  very 
urgent  and  interrelated  problems.  The  near-ground  sound  propagation  has  its  peculiarities  and  differs  from  the 
sound  propagation  in  free  space.  In  addition  to  the  meteorological  parameters  (the  air  temperature,  the 
atmospheric  pressure,  the  air  humidity,  and  the  wind  speed  and  direction),  geometric  factors,  for  example,  mutual 
arrangement  of  a  source,  a  receiver,  and  an  underlying  surface  also  affect  the  near-ground  sound  propagation.  In 
some  cases,  these  factors  may  dominate  in  comparison  with  the  effect  of  the  meteorological  parameters  and 
hence  call  for  much  attention.  In  the  present  project,  the  above-indicated  problems  have  been  studied  both 
theoretically  and  experimentally.  The  results  of  investigations  were  presented  in  7  interim  reports  on  this 
projects.  This  final  report  presents  the  results  of  our  investigations  into  the  effect  of  the  atmospheric  channel  on 
the  acoustic  wave  propagation  above  the  ground  as  well  as  of  investigations  into  the  amplitude  and  phase 
fluctuations  of  acoustic  waves  propagating  through  the  turbulent  atmosphere  near  the  ground  and  the  structure  of 
the  temperature  and  wind  velocity  fields  and  the  vertical  profiles  of  the  outer  scale  of  atmospheric  turbulence. 


P 


* 
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2  INFLUENCE  OF  THE  ATMOSPHERIC  CHANNEL  ON  THE  SOUND 
PROPAGATION  ABOVE  THE  GROUND  SURFACE 


It  is  well  known  that  the  atmospheric  sound  recorded  from  a  stationary  source  essentially  depends  on 
meteorological  conditions.  This  is  caused  by  the  high  sensitivity  of  the  total  energetic  losses  of  sound 
to  such  meteorological  parameters  as  wind,  temperature,  relative  air  humidity,  and  atmospheric 
pressure  Mean  values  of  humidity,  temperature,  and  pressure  determine  not  only  the  coefficient  of 
molecular  absorption  of  the  acoustic  radiation  at  a  fixed  frequency,  but  also  its  frequency  dependence 
(Brown  E  H  and  Hall  F.F.,  1978;  Harris  C.M.,  1966).  Fluctuations  of  acoustic  refractive  index,  caused 
by  turbulent  pulsations  of  the  meteorological  parameters,  affect  an  interference  pattern  of  directly 
transmitted  and  reflected  acoustic  rays  and  cause  the  excess  turbulent  attenuation  (Brown  E.H.  and 
Hall  F  F  1978-  Ingard  U.,  1953).  In  case  of  near-ground  sound  propagation  at  distances  more  than 
1  km,  sound  attenuation  is  primarily  affected  by  refraction  on  gradients  of  temperature  and  wind 
(Ostashev  V.E.,  1992),  resulting  in  the  waveguide  or  nonwaveguide  regime  of  sound  propagation.  In 
the  first  case,  rays  bend  downward  and  hence  undergo  multiple  reflections  from  the  ground.  Relatively 
small  values  of  sound  attenuation  are  typical  of  this  regime  of  sound  propagation.  In  the  second  case, 
rays  bend  upward,  and  a  shadow  zone  is  formed  near  the  ground  at  certain  distances  from  the  source. 
Only  very  weak  sound,  caused  by  turbulent  scattering  in  the  upper  atmospheric  layers,  penetrates  into 
this  zone. 


Previously  we  described  the  application  package  Outdoor  Acoustics  (Abramov  N.G.  et  al.,  1994; 
Bogushevich  A.Ya.  and  Krasnenko  N.P.,  1993,  1996)  for  calculation  of  the  sound  pressure  level  in 
case  of  sound  propagation  through  the  atmosphere  in  the  waveguide  or  nonwaveguide  regime  at 
distances  up  to  10  km  and  reported  the  field  measurements  that  confirmed  reliability  and  accuracy  of 
the  calculated  data.  In  the  present  work,  an  attempt  is  made  to  analyze  numerically,  with  the  use  of 
this  application  package,  the  effect  of  the  meteorological  parameters  on  audibility  of  sound  in  the 
atmosphere.  The  most  interesting  results  of  calculations  are  shown  in  Figures  1-4. 


The  dependence  of  the  total  attenuation  L  on  the  wind  component  parallel  to  the  direction  of  sound 
propagation  \4>  is  shown  in  Fig.  1.  Here,  the  waveguide  regime  of  propagation  is  realized  for  Vp>0  and 
the  nonwaveguide  regime  -  for  l/p<0.  Refraction  by  wind  velocity  gradients  always  results  in  the 
increase  of  sound  attenuation.  In  the  waveguide  regime  of  propagation,  the  increase  of  Vp  is 
accompanied  by  a  relatively  small  increase  in  |/|,  as  a  rule,  by  no  more  than  6-10  dB,  connected  with 
the  change  of  the  beam  focusing  pattern  compared  with  the  conventional  spherical  beam  divergence. 
At  the  same  time,  in  the  nonwaveguide  regime  of  propagation,  the  excess  sound  attenuation  due  to 
wind  is  much  greater  and  may  even  exceed  100  dB.  In  this  case,  the  sharp  increase  of  |Z|  by  15- 

20  dB  is  first  observed  as  j  V^|  increases  and  the  observation  point  enters  the  acoustic  shadow  zone. 


With  the  further  increase  of 


the  sound  attenuation  continues  to  increase  noticeably,  because  the 


atmospheric  layer,  from  which  a  scattered  signal  comes  to  the  observation  point,  is  displaced  upward. 
The  change  of  the  distance  from  the  source  d  affects  insignificantly  the  form  of  the  dependence  L(VP), 
but  shifts  noticeably  the  absolute  values  of  L 
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Fig.  1.  The  dependence  of  sound  attenuation  in  the  atmosphere,  in  dB,  on  the  wind  component  Vp 
parallel  to  the  direction  of  sound  propagation  for  P-750  mm  Hg,  U-75%,  T=20°C,  y=0 /km, 
C2-10'6 m213,  MOOO  Hz,  and  distances  to  the  source  d=1(1),  3(2),  6(3),  and  10  km(4). 
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Fig.  2.  Sound  attenuation  in  the  atmosphere,  in  dB,  as  a  function  of  the  wind  direction  F  for  U=75%, 
T=20°C,  g=0°/km,  Cn2=10'6 rnm,  f=1000 Hz,  and  V=\v\  = 0 (1),  5  (2),  10  (3),  and  15 m/s  (4). 


L,  dB 


Figure  2  shows  sound  audibility  as  a  function  of  the  wind  direction  F  counted  off  from  the  source- 
observer  direction.  Curve  1  was  calculated  without  refraction  (V=0).  It  can  be  easily  seen  that  at 
V=const  the  change  of  F  is  equivalent  to  the  change  of  the  wind  component  Vp  parallel  to  the  direction 
of  sound  propagation.  Therefore,  curves  2,  3,  and  4  exhibit  practically  the  same  behavior  as  in  Fig.  1. 
From  Fig.  2  it  follows  that  the  wind  component  perpendicular  to  the  direction  of  sound  propagation 
affects  insignificantly  its  audibility. 

Figure  3  shows  the  dependence  of  L  on  Vp  for  indicated  values  of  the  structure  parameter  of  refractive 
index  fluctuations  Cn2that  characterize  the  intensity  of  the  atmospheric  turbulence.  As  expected,  Cn 
affects  most  strongly  the  sound  pressure  level  in  the  nonwaveguide  regime  of  sound  propagation.  In 
this  case,  the  recorded  pressure  level  is  almost  proportional  to  Cn2.  Because  Cn  varies  over  a  wide 
range,  sound  audibility  also  undergoes  large  variations.  In  our  case  experimental  values  changed  from 
2  to  5  km. 


8 


L,  dB 


Fia  3  Sound  attenuation  in  the  atmosphere,  in  dB,  as  a  function  of  the  wind  component  Vp  parallel  to 
the  direction  of  sound  propagation  for  U=75%,  T=20°C,  g=0°/km,  f=1000Hz,  d=6  km,  and 
Cn2=10'a  (1),  10'7  (2),  10 f*  (3),  and  10's  m2*  (4). 


In  the  waveguige  regime  of  propagation,  the  observed  sound  pressure  level  is  a  sum  of  pressure  levels 
of  individual  rays  propagating  along  paths  of  different  lengths.  The  total  pressure  level  differs 
qualitatively  as  a  function  of  amplitude  and  phase  fluctuations  of  rays.  For  small  values  of  C„ ,  the 
amplitude  and  phase  fluctuations  of  rays  are  also  small  and  hence  can  be  partly  coherent.  In  this 
case,  because  the  ray  path  lengths  depend  on  the  V  profile,  the  linear  increase  of  V?  results  in  the 
interference  sound  pressure  pattern  vividly  seen  in  Fig.  3  for  curves  1  and  2.  For  large  values  of  Cn , 
the  rays  are  incoherent  and  interference  is  not  observed.  Coherence  starts  to  break  at  smaller  values 
of  Cn2  as  the  distance  between  the  source  and  the  receiver  increases. 


Fig  4  Sound  attenuation  in  the  atmosphere,  in  dB,  as  a  function  of  the  temperature  gradient  g  for 
U-75%,  T=2CfC,  C2=10 r*  m2/3,  M000  Hz,  and d=1  (1),  3  (2),  6  (3),  and  10 km  (4). 

The  effect  of  temperature  gradient  on  sound  attenuation  in  the  atmosphere  is  illustrated  by  Fig.  4.  For 
g>0  the  waveguide  regime  of  sound  propagation  is  realized,  as  in  the  case  of  wind  refraction  for  Vp>0, 
but  in  contrast  to  it,  additional  beam  focusing  takes  place  resulting  in  sound  amplification.  Therefore, 
the  presence  of  temperature  inversion  in  the  atmosphere  increases  the  sound  level.  In  calculations, 
the  sound  source  was  at  the  distance  r,=  5  m  above  the  ground  and  the  observation  point  was  at 
r2=1.5  m.  In  the  nonwaveguide  regime  of  sound  propagation,  the  zone  of  direct  sound  propagation 
always  exists  in  the  vicinity  of  the  source  followed  by  the  acoustic  shadow  zone.  For  hight  source  and 
receive  the  distance  to  this  shadow  zone  is  determined  by  the  curvature  radii  of  ray  paths,  which  are 
noticeably  smaller  in  case  of  refraction  by  temperature  gradients  than  in  case  of  refraction  by  wind 
gradients.  Therefore,  when  g  <0  and  the  observation  point  is  in  the  shadow  zone,  the  sharp  increase 
of  sound  attenuation  takes  place  for  large  absolute  values  of  g.  The  further  increase  of  the  sound 
attenuation  in  the  shadow  zone  as  |g|  increases  is  caused  by  the  upward  displacement  of  the  lower 
boundary  of  the  region  of  sound  scattering,  as  in  the  case  of  wind  refraction. 
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3.  ANISOTROPY  OF  THE  AMPLITUDE  FLUCTUATIONS  OF  THE 
ACOUSTIC  WAVES  PROPAGATED  THROUGH  THE  TURBULENT 
ATMOSPHERE  NEAR  GROUND 


The  amplitude  fluctuations  of  an  acoustic  wave,  propagating  in  the  ground  layer  of  the  turbulent 
atmosphere  under  conditions  of  reflection  from  the  Earth’s  surface,  were  studied  theoretically 
(Diaafile  CA  1980-  Clifford  S.F.,  Lataitis  R.L.,  1983;  Hidaka  T.,  Kageyama  K.,  Masuda  S., 
1985)Sand  experimentally  (Parkin  P.H.,  Scholes  W.E.,  1965;  Daigle  G. A.,  Piercy  J.E.,  Embleton 
T  F.W.,  1978;  Bochkarev  N.N.,  Krasnenko  N.P.,  1986).  In  theory  the  peculiarities  of  an  acoustic 
wave  propagation  in  a  randomly  inhomogeneous  moving  medium  (Ostachev  V  E.  1992)  were  no 
taken  into  account  while  experimental  investigations  were  conducted,  as  a  rule,  at  a  fixed  path 
orientation  relative  to  the  wind  velocity  vector.  In  the  experiment  conducted  we  have  developed 
and  refined  investigations,  described  in  (Patrushev  G.Ya.,  Postov  A.P.,  1996). 

The  difference  of  the  results  presented  in  this  report  from  those  published  in  (Patrushev  G.Ya 
Postov  A  P  1996)  is  as  follows.  We  used  more  precise  equipment.  In  the  experiment,  descn  e 
in  (Patrushev  G.Ya.,  Postov  A.P.,  1996),  for  different  directions  different  types  of  microphones 
and  loudspeakers  were  used.  There  was  also  used  an  analog  detecting  of  signals  with  their 
following  smoothing  by  low-frequency  analog  filters  of  the  second  order  with  further  discretization 
at  a  frequency  of  100  Hz.  This  could  introduce  additional  error  in  the  estimation  of  spatial 
anisotropy  of  the  acoustic  wave  propagation.  In  this  study  we  have  used  quite  identical  measuring 
condenser  microphones  and  dynamic  speakers.  A  2-channel  12-bit  high-speed  ADCs  (conversio, 
time  was  3  ms)  were  used  for  digitizing  signals  and  correspondingly  very  high  frequency  of 
discretization  of  the  signals  themselves  from  the  microphones.  Also  we  used  practically  new  method 
of  the  amplitude  separation.  The  fast  computer  enabled  using  the  programs  working  m  real  time 
and  giving  more  exact  results.  All  this  allowed  us  to  reach  the  identity  of  the  through  channels 
(emitter-microphone-recording  accurate  to  1%  in  the  operating  frequency  band). 

The  approach  to  experimental  data  processing  using  low-frequency  filtration  with  the  following 
decimation  (without  smoothing)  is  quite  different  than  that  used  in  (Patrushev  G.Ya.,  Postov  A.P., 
1996)  and  it  allowed  more  exact  study  of  the  low-frequency  portion  of  the  spectrum. 

In  contrast  to  the  study  presented  in  (Patrushev  G.Ya.,  Postov  A.P.,  1996)  we  arranged  full 
meteorological  measurements. 

We  analyzed  the  data  acquired  at  different  wind-temperature  stratification  whereas  in  (Patrushev 
G  Ya  Postov  A  P  1996),  the  data  obtained  only  at  one  (unstable)  stratification  were  available. 
The  synchronous  measurements  we  have  conducted  of  the  relative  amplitude  fluctuations  of  an 
acoustic  wave  show  significant  anisotropy  of  fluctuations,  depending  on  the  average  wind  velocity, 
other  conditions  being  the  same.  The  diagram  of  the  experiment  is  shown  in  Fig.  5. 


FIG.  5.  Block-diagram  of  the  experiment. 
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Two  quasi-spherical  waves  directed  along  the  average  wind  velocity  and  perpendicularly  to  it  were 
measured  simultaneously.  Receiving  microphones  were  at  L  ~  30  m  from  the  emitting  horns. 
Emitters  and  receivers  were  located  at  the  height  h-  1.5  m.  The  loudspeaking  horns  with  an 
output  of  50  W  used  as  the  emitters,  and  as  the  receivers  were  used  measuring  condenser 
microphones  MV201  by  RFT,  with  the  antiwind  protection  shields.  Each  microphone  was 
connected  to  an  ac  pre-amplifier.  The  width  of  the  receiver's  directional  pattern  was  close  to  GO", 
and  receiving  microphones  have  the  directional  pattern  of  the  cardioid  shape.  Signals  from  the 
amplifiers  of  the  microphones  entered  the  third-octave  filters  with  the  central  frequency  at 
1.25  kHz.  The  output  signals  from  the  filters  were  digitized  with  a  12-bit  ADC  at  a  sampling 
frequency  of  160  kHz  and  after  the  amplitude  detection  their  values  were  entered  in  the  computer 
memory  at  a  frequency  of  100  Hz.  The  amplitude  detection  was  organized  programmly  using  the 
following  algorithm:  1600  digital  values  of  the  ADC  codes  were  memorized  in  a  buffer,  then  they 
were  detected  and  the  maximum  value  which  was  recorded  as  the  instantaneous  value  of  the 
amplitude  was  found  from  the  obtained  data  array.  This  algorithm  was  realized  on  the  assembler 
using  the  register  addressation,  thus  enabling  its  real  time  operation. 

Instantaneous  values  of  the  wind  and  temperature  components  were  measured  with  a  three- 
component  ultrasonic  anemometer-thermometer  located  at  the  height  of  h  ~  1.5  m  at  the  bisectrix  of 
the  angle  between  the  microphones.  The  flat  field  of  flattened  oats  15—20  cm  high  was  the 
underlying  surface.  The  frequency  of  the  emitted  wave  was  f  -  1.25  kHz.  Measurements  were 
conducted  both  at  stable  and  unstable  stratification  of  the  atmosphere  and  the  time  of  an  individual 
sampling  was  300  s.  The  range  of  variation  of  the  average  temperature  during  the  measurements  at 
the  height  of  h  ~  1.5  m  was  1.5-6. 8°,  relative  humidity  of  the  air  99-75%,  and  the  pressure  752- 
755  mm  of  the  mercury  column. 

Further  statistical  processing  of  the  obtained  time  series  of  the  instantaneous  values  of  the 
amplitude  included  the  removal  of  rough  errors  using  Chcbyshev  inequality,  low  frequency 
filtration  using  Butterwort  filter  of  the  6th  order  (Otrtes  R.,  Enokson  L.  1982)  and  agreed  with  it 
decimation  of  the  data  with  the  index  8.  Then  the  removal  of  the  trend  was  concluded  and 
autospectra  were  calculated  after  the  FFT  processing.  The  effect  of  infiltration  was  removed  by  the 
GEO  window  (Otnes  R.,  Enokson  L.  1982)  after  the  calculation  of  the  Fourier  coefliciont.  In 
parallel  with  the  autospectra  the  average  values  of  the  amplitudes,  their  variances  and 
autocorrelation  functions  were  also  calculated. 

„,_(A2)-(A)2  (a2) 

P  2  ~  2 

(A)  (A) 


’  0  1  2  3  4 


FIG .  6.  The  dependence  of  the  relative  rms  deviation  p  on  the  average  wind  velocity  V  for  the  N 

paths  along  (o)  and  across  (A)  the  wind. 

In  Fig.  6  are  shown  the  values  of  the  relative  rms  deviation  p  of  the  amplitude  fluctuations  A, 
depending  on  the  value  of  the  average  wind  velocity  V  for  longitudinal  and  transverse  relative  to 
the  wind  direction  paths.  The  straight  lines,  constructed  as  the  best  linear  approximation  ot  the 
experimental  data  are  also  shown  in  the  figure.  As  the  obtained  data  show  much  stronger 
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dependence  of  the  amplitude  fluctuations  on  the  wind  speed  is  observed  for  the  cross  wind  path  as 
compared  to  lliat  lor  the  path  along  the  wind. 


b 

I'lG.  7.  The  temporal  spectra  ( a)  and  autocorrelation  functions  (b)  oj  the  amplitude  fluctuations 
oj  an  acoustic  wave  for  the  paths  atony  ( 1 )  and  across  (2)  the  wind  at  a  stable  stratification  of 
the  atmosphere  (£,  -  +0.1):  1  -  Vi  =  0  008  m/s;  ax  =  0.3  n»/s;  2  -  Vx  =  0.9  m/s;  a,  =  0.4  m/s. 


FIG.  <S\  The  temporal  spectra  (a)  and  autocorrelation  functions  (b)  of  the  amplitude  fluctuations 
of  an  acoustic  wave  for  the  paths  along  ( 1)  and  across  (2)  the  wind  at  an  unstable  stratification 
of  the  atmosphere  <$  -  -0.9):  /  -  Vx  =  0.0(i  m/s;  aa  =  0.27  m/s;  2  -  VL  =  1.2  m/s; 

—  0.32  m / s. 

In  hies.  7  and  8  are  shown  the  dimensionless  temporal  spectra  of  the  amplitude  fluctuations  for 
both  paths  and  corresponding  to  this  spectra  normalized  temporal  autocorrelation  functions 
obtained  at  stable  (I'ig.  7)  and  unstable  (Fig.  8)  stratifications  of  the  atmosphere.  I  he  scale  o 
similarity  considering  variations  of  the  average,  perpendicular  to  the  wind  velocity  Vu  and 
fliictuational  component  oL,  is  the  characteristic  frequency  fQ=  |(l/2L  +  oi )/ (X./.)|,/2  (Patrushev 
G.Ya.,  Post ov  A.P.,  1994),  where  X  is  the  length  of  the  wave  of  radiation,  I.  is  the  length  of  the 
path.  From  the  data  presented  one  can  see  that  the  time  spectrum  for  the  path  perpendicular  to  the 
average  wind  undergoes  broadening  to  the  low-frequency  region  as  compared  to  the  tunc  spectrum 
obtained  for  the  path,  directed  along  the  average  velocity  ol  the  wind,  and  moreover,  under 
unstable  stratification  the  broadening  is  considerably  more  noticeable  and  forms  approximately  hall 
an  order  of  magnitude  on  the  scale  of  the  relative  frequences. 
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4.  PHASE  FLUCTUATIONS  OF  THE  ACOUSTIC  WAVE 
PROPAGATING  NEAR  GROUND 


Fluctuations  of  phase  difference  of  acoustic  waves  propagating  in  the  turbulent  atmosphere  were  theoretically 
studied  in  detail  in  (Tatarskii  V.  L,  1961,  Chernov  L.A.,  1975),  where  their  linear  dependence  on  the  structural 
characteristic  of  the  refractive  index  fluctuations  Q2  and  on  the  length  of  the  propagation  path  L  and  I716 
dependence  on  the  wave  number  were  established.  The  temporal  spectrum  of  phase  difference  fluctuations 
depends  on  the  length  of  measurement  path  and  the  average  wind  velocity  component  perpendicular  to  the 
propagation  path.  The  form  of  the  probability  density  of  phase  difference  fluctuations  /(<t>)  was  not  discussed; 
however,  in  the  geometric  acoustics  approximation  and  in  the  context  of  the  smooth  perturbation  method  the 
Gaussian  probability  density  of  phase  difference  fluctuations  should  be  expected  when  the  central  limit  theorem 
is  applicable.  Up  to  now,  the  problem  of  saturation  of  phase  fluctuations  under  different  propagation  conditions 
also  was  not  discussed. 

Experiments  performed  under  conditions  when  the  effect  of  the  ground  can  be  neglected  agree  fairly  well  with 
the  theory  based  on  the  method  of  geometric  acoustics  for  the  Kolmogorov  spectrum  of  atmospheric  turbulence. 
Experimental  data  on  the  phase  difference  fluctuations  a*2  and  on  the  fluctuations  of  the  mutual  coherence 
function  were  reported  in  (Daigle  G.A.  et  al.,  1983).  The  measurements  there  were  carried  out  for  propagation 
paths  of  lengths  L  from  2  to  300  m  and  were  accompanied  by  independent  measurements  of  the  atmospheric 
meteorological  parameters.  Two  transient  phase  meters  were  used  to  measure  the  phase  difference.  They  were 
capable  of  measuring  phase  differences  up  to  2n  rad  and  were  shifted  by  n  rad.  It  should  be  noted  that  operation 
of  the  phase  meters  was  controlled  by  a  computer,  which  assigned  the  constant  phase  difference  to  those  readings 
of  the  phase  meter  for  which  the  input  signal  died  down  and  was  less  than  the  threshold  specified  by  the  operating 
instructions  to  these  phase  meters.  Obviously,  this  can  be  considered  as  a  disadvantage,  because  the  signal  phase 
undergoes  fast  variations  exactly  when  the  input  signal  dies  down.  These  fast  phase  variations  cannot  be 
measured  by  this  method.  Nevertheless,  Daigle  et  al.,  1983  obtained  good  agreement  between  their  experimental 
data  and  calculations  for  the  Gaussian  spectrum  of  atmospheric  turbulence  with  the  finite  value  of  the  outer  scale 
of  turbulence  Lo.  It  should  be  mentioned  that  the  measurements  were  carried  out  above  asphalt-coated  underlying 
surface,  and  the  receiver  was  located  outside  the  refraction  shadow  zone.  To  this  end,  the  receiver  height  above 
the  ground  was  changed  with  the  increase  of  the  measurement  path  length.  The  increase  of  the  receiver  height 
should  cause  the  increase  of  the  outer  scale  of  atmospheric  turbulence;  however,  for  their  calculations  Daigle 
et  al.,  1983  used  constant  value  of  the  outer  scale  of  turbulence  Zo=l  .1  m  even  for  the  receiver  height  as  great  as 
7  m.  Other  characteristics  of  the  acoustic  wave  phase  fluctuations  were  not  considered  in  that  paper. 

We  measured  the  acoustic  wave  phase  difference  for  propagation  paths  50,  100,  200,  and  400  m  long.  Receiving 
microphones  were  spaced  at  tfel  m.  Their  heights  above  the  ground  were  /?i~l  in  in  the  first  series  ol 
measurements  and  hf* 0  in  the  second  series  of  measurements,  that  is,  the  receivers  were  placed  directly  on  the 
ground.  A  source  was  at  heights  4=(1.3~1.5)m  above  the  ground.  For  this  geometry  of  the  experiment,  the 
source  can  be  considered  as  a  source  of  spherical  acoustic  waves  with  a  frequency  of  1  kHz. 

Figure  9  shows  the  rms  values  of  the  phase  difference  fluctuations  versus  the  parameter 
p=0.308  C^l)  ,/6(27tJ?  Vs)116  that  specifies  the  conditions  of  sound  wave  propagation  in  free  space  as  functions  of 
the  structural  characteristic  of  the  refractive  index  fluctuations  propagation  path  length  £,  and  the  acoustic 
wave  number  Here,  Vs  is  the  speed  of  sound.  In  the  figure,  circles  are  for  the  first  series  of 

measurements  with  A=1  m,  and  crosses  are  for  the  second  series  of  measurements,  when  the  receivers  were 
placed  on  the  ground.  During  our  measurements  we  could  not  record  signals  simultaneously  for  different 
propagation  path  lengths  L  under  the  same  meteorological  conditions  (determined  by  the  structural  characteristic 
of  the  refractive  index  fluctuations  and  gradients  of  the  air  temperature  and  wind  velocity).  Therefore,  it  was 
quite  reasonable  to  choose  the  parameter  p  to  characterize  the  conditions  of  acoustic  wave  propagation.  It  can  be 
seen  from  the  figure  that  the  rms  values  of  phase  difference  increase  when  the  parameter  P  increases  by  about 
two  orders  of  magnitude.  On  the  semi-logarithmic  scale,  this  dependence  can  be  approximated  by  a  straight  line. 
However,  some  data  recorded  for  small  p  disobey  this  dependence.  It  should  be  noted  that  for  long  propagation 
paths  with  £=400  m  the  experimental  phase  difference  a*  was  close  to  that  reported  in  Daigle  et  al,  1983.  We 
could  not  compare  our  data  more  accurately,  because  (Daigle  et  al,  1983)  did  not  indicate  the  distance  between 
the  receivers.  It  also  can  be  seen  from  the  figure  that  the  behavior  of  the  phase  difference  is  practically  the  same 
for  the  receivers  located  at  two  heights  h=  0  and  1  m. 
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The  probability  density  characterizes  more  adequately  the  fluctuations  of  a  random  variable.  IJp  to  now,  as 
already  mentioned  above,  no  experimental  data  were  reported  for  this  characteristic  (at  least,  they  arc  unknown  to 
us).  This  characteristic  also  was  not  investigated  theoretically.  At  the  same  time,  our  experimental  procedure  can 
be  used  to  study  the  probability  density  of  the  phase  difference  /3(<|)).  Curve  I  in  big.  10  shows  the  histogram  of 
the  phase  difference  /{ <|>)  for  the  propagation  path  of  length  £-100  m,  height  of  the  receivers  h  \  m,  g^-0.44, 
and  [5-  0.037.  Curve  2  shows  the  Gaussian  probability  density  drawn  for  the  same  variance.  It  can  be  seen  from 
the  figure  that  the  distribution  of  the  measured  fluctuations  of  the  acoustic  wave  phase  difference  significantly 
deviates  from  the  normal  distribution,  namely,  its  modal  value  is  more  probable  and  large  values  of  phase 
difference  are  less  probable.  Such  behavior  of  the  probability  density  was  typical  of  practically  all  experimental 
histograms  of  phase  difference  fluctuations,  not  only  for  measurements  with  the  receivers  located  directly  on  the 
ground,  but  also  for  measurements  with  the  receivers  placed  at  heights  of  1  m  above  the  ground  and  spaced  in  the 
horizontal  direction. 
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Fig.  10.  Comparison  of  the  experimental  probability  density  /?<(>)  (curve  1) 
with  the  norma!  probability  density  (curve  2) 


Figure  11  shows  the  temporal  spectra  of  fluctuations  of  acoustic  wave  phase  difference  recorded  with  Ihe 
receivers  located  at  heights  of  1  m  above  the  ground  under  conditions  of  strongest  wind  velocity  VtA  m/s 
measured  at  a  height  of  3  m  above  the  ground  (curve  2)  and  when  the  wind  velocity  measured  at  a  height  of  3  m 
was  ^=0.39  m/s  (curve  1).  The  logarithm  of  the  frequency  /(in  Hz)  is  plotted  on  the  abscissa,  and  the  spectral 
power  density  G( $  normalized  to  the  variance  a*2  and  multiplied  by  the  frequency  /is  plotted  on  the  ordinate. 
The  dimensionless  spectrum  ct$2  can  be  used  to  estimate  the  contribution  of  different  frequencies  to 

the  fluctuation  variance  in  a  wide  spectral  range.  As  follows  from  the  theory  of  acoustic  wave  propagation,  the 
temporal  spectrum  of  phase  fluctuations  broadens  linearly  with  the  increase  of  the  wind  velocity  component 
orthogonal  to  the  propagation  path  in  the  free  atmosphere  (the  similarity  parameter  is  fd!V l).  As  can  be  seen 
from  the  figure,  a  10-fold  increase  of  the  wind  velocity  component  perpendicular  to  the  propagation  path  caused 
the  shift  of  the  position  of  the  spectral  power  density  maximum  toward  larger  frequencies.  Its  frequency  increased 
no  more  than  twice.  That  is,  the  dependence  of  the  temporal  spectrum  of  the  phase  difference  fluctuations  on  the 
perpendicular  component  of  the  wind  velocity  is  very  weak  in  case  of  near-ground  propagation.  It  also  should  be 
noted  that  the  characteristic  fluctuation  frequencies  were  no  more  than  20-30  Hz,  and  the  maximum  of  the 
dimensionless  spectrum  was  about  0.3,  which  was  typical  of  the  propagation  in  the  free  atmosphere.  An  analysis 
of  the  temporal  phase  difference  spectra  showed  that  they  did  not  depend  on  the  propagation  path  length  and  the 
intensity  of  atmospheric  turbulence  for  the  considered  heights  of  receivers  (/r=0  and  1  m)  when  all  other  i actors 
remained  unchanged. 
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To  calculate  the  mutual  coherence  function,  v/e  used  the  discrete  Hilbert  transform.  According  to  the  definition, 
let  Aj(/)  be  the  instantaneous  value  ol  the  signal  \|(/)  .4i(/)exp(i(27i//iai(0)|,  where  A\(f)  and  cxi(/)  are  leal 
random  functions,  and  analogously,  a,</)  /I  >(/)exp[/(2n/?i  (x>(/))|.  Then  the  mutual  coherence  function  Ci,2  can  be 
written  as 

^ _ <A\(0a\(0> _ 

[<  Ai(f)  A\ (0  >H[<  Aid)  AiiO  >\~ 

where  (he  angular  brackets  denote  time  averaging  defined  by  the  formula 

(j 

Figure  12  shows  tire  mutual  coherence  functions  on  the  propagation  path  of  length  50  m  for  the  values  of  the 
parameter  P“0 .009-0.01  and  for  the  propagation  path  of  length  400  m  and  (1  0  9-1.6.  It  can  be  seen  that  the 
degree  of  coherence  decreases  noticeably  for  longer  propagation  path.  It  also  should  be  noted  that  the  mutual 
coherence  function  measured  in  (Daigle  et  al,  1983)  at  a  frequency  of  I  kllz  for  the  propagation  path  of  length 
100  ,n  with  the  receivers  spaced  at  a  distance  of  I  m  was  much  greater.  This  may  be  explained  by  the  tact  that  the 
receivers  were  placed  at  greater  heights  above  the  ground  than  in  our  measurements  as  well  as  by  different  type 
ofthc  underlying  surface  (hard  covering  and  cultivated  field  sowed  with  peas  and  oats). 
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Fig.  12.  Mutual  coherence  function  Tu  on  the  paths  50  (P-0-01)  and  400  m  long  (p=0.9-1.6) 


To  summarize  the  results  of  our  investigations  of  statistical  characteristics,  we  can  conclude  the  following: 

In  summer  the  nns  values  of  the  intensity  fluctuations  of  acoustic  waves  propagating  near  the  ground  varied  by 
[bout  an  order  of  magnitude  on  paths  of  lengths  from  50  to  400  m  for  the  receivers  spaced  at  1  m  as  functions 
the  intensity  of  atmospheric  turbulence. 

The  probability  density  of  fluctuations  of  phase  difference  significantly  deviated  from  the  normal  one.  The 
temporal  spectra  of  fluctuations  depended  on  the  wind  velocity  measured  at  the  height  of  receiver  location  much 
weaker  thm  in  case  of  line-of-sight  propagation.  With  the  increase  of  the  propagation  path  length  the  mutual 
coherence  function  measured  at  a  frequency  of  1  kHz  decreased.  For  path  400  m  long  and  receivers  spaced 
m  it  was  between  0-0.3. 

Our  measuremem  method  is  rather  adaptable.  It  allowed  us  to  study  all  statistics  of  phase  fluctuations.  To 
improve  the  method,  we  should  improve  the  measurement  procedure  including  synchronous  recording  of signa is 
for5 several  propagation  paths  and  different  separations  of  the  receivers  and  measurements  of  the  outer  scale  of 

atmospheric  turbulence. 
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5.  EXPERIMENTAL  VERIFICATION  OF  A  MODEL  FOR  THE  SPECTRAL 

tensor  of  the  wind  velocity  field 


Detailed  structure  of  small-scale  turbulence  of  lower  layers  of  the  atmosphere  is  interesting  from 
Sent  pointsof  view.  There  is  a  wide  scope  of  applied  problems  requiring  adequate  descnpt, on 
of  the  spatial  spectrum  of  turbulent  pulsations  of  wind  velocity  in  the  surface  layer  nformation  of 
such  a  kind  is  necessary,  for  instance,  in  the  problems  on  propagation  of  optical  and  radio  waves  in 
the  atmosphere  in  solving  problems  of  acoustic  and  lidar  sounding,  for  taking  account  of  pollutant 
ecological  problems,  for  estimation  of  spatial  variability  of  wind  pressure  in  calculating 

constructional  stresses,  etc. 

According  to  modern  experimental  and  theoretical  data,  the  three-dimensional  spatial  spectrum  of 
turbulent  pulsations  of  wind  velocity  is  a  tensor  of  the  second  rank. 

In  the  approximation  of  isotropic  turbulence  (Tatarskii  V.I.,  1961),  the  spectral  tensor  Oi;(k)  is 
completely  determined  by  a  single  scalar  function  of  the  wave  number  k  -  Ik  I 

<Di;-(k)  =  <4icfe2)-«{5iy  -  k,  kj  ,  (51> 

namely,  by  the  energy  spectrum  »(k)  for  which  the  Kolmogorov- Obukhov  law  of  •‘five  ihi.ds”: 
E(k)  ~  e2/3/r5  /;',  where  £  is  the  dissipation  rate  of  the  turbulence  kinetic  energy,  should  hold  fo 
reasons  of  dimensionality. 

The  isotropic,  approximation  in  the  surface  layer  sufficiently  well  works  in  the  inertial  interval  of 
the  wave  number.  However,  l he  turbulence  becomes  anisotropic,  with  the  increase  of  '‘ales 

and  approaching  the  energy  interval  of  the  spectrum  (Newstadl  I  .I.M.,  Van  Dop  H.,  t.)8a).  i 
fact  when  the  distance  between  the  observation  points  of  wind  velocity  fluctuations  is  comparable 
S  the  outer  turbulence  scale,  the  hypothesis  on  the  spatial  isotropy  of  fluctuations  tads.  So  it  » 
necessary  to  construct  better  models  for  the  spectral  tensor. 

According  to  the  Monin-Obukhov  similarity  theory,  the  dimensionless  spectra  of  fluctuations  of 
wind  velocity  vector  components  in  the  surface  layer  of  the  atmosphere  at  the  height  *  are  universal 
functions  of  the  dimensionless  wave  number  kz  and  the  stratification  parameter  z/L  where  L  is  the 
Obukhov  scale.  This  means  that  specified  characteristic  frequencies  corresponding  to  spectra 
maxima  low-frequency  boundaries  of  the  inertial  interval,  etc.  depend  on  the  stratification 
(Zilitinkevich  S.S.,  1970).  Thus,  the  similarity  theory  permits  one  to  make  some  concrete 
conclusions  about  the  shape  of  the  spectra  under  different  states  of  the  surface  layer  stability. 
Thus,  the  attempts  to  take  into  account  some  universal  regularities  m  the  surface  layer  by  the 
similarity  theory  and  experimental  data  are  quite  urgent  in  simulation  of  the  spectral  tensor. 

The  kinematic  model  of  turbulence  developed  by  Kristensen  L„  Lenshow  D.H.,  Kirkegaard  P 
Courtney  M„  1989  is  of  special  interest  among  the  investigations  on  this  problem.  It  is  based  on 
tire  representation  of  the  spatial  spectrum  of  the  velocity  field  fluctuations  as  an  an.so  rop.c  tensor 
Or(k).  The  following  form  is  proposed  in  the  above-mentioned  paper  for  the  spectral  tensor  in  e 
approximation  of  a  homogeneous  incompressible  turbulent  flow: 


d)y(k)  =  Yj  A<*>  {8/i  -  k,  hi  r2}{8,,  -  k,  kj  r2}. 

1=1 


(5.2) 


In  contrast  to  the  case  of  isotropic  turbulence  (5.1),  the  description  of  the  tensor  (5.2)  needs  three 
independent  real  scalar  functions  /!,(*).  A2(k),  As(k)  determining  the  energy  of  turbulent  vortices 
along  three  orthogonal  directions  assigned  by  unit  vectors  it,  »2.  '3- 


Variances  al ,  al,  al,  integral  scales  /„,  lv,  lw,  and  dimensionless  parameters  pu,  P,,  ^ 
characterizing  the  spectrum  inclination  in  the  transition  range  from  the  inertial  interval  to  the 
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energy  interval  of  wave  numbers  are  used  in  (Kristensen  L.,  Lenshow  1L1L,  Kirkegaard  P., 
Courtney  M.f  1989)  as  the  parameters  describing  the  behavior  of  the  spectra  of  longitudinal, 
transverse,  and  vertical  components  u,  w  of  the  velocity.  Here,  the  analytical  expressions  for 
one  dimensional  spatial  spectra  of  fluctuations  of  the  velocity  components  have  the  form 


/'«(*)  = 


_  luCTu 


i  +  i_4A_ 


2f‘u 
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(,f‘u 


l-vu<) 


lo  <Jo 
2n 
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(5.3) 


where  a(p)  =  {jr^tTlS /6j.i |}/{l'(  1 /2ja)r[ ( 1  /3fi) ]}  is  a  dimensionless  constant  determined  by  the 

oO 

normalization  condition  a2  =  J  F(k)dk  for  each  concrete  value;  k  is  the  wave  number  along  the 

—co 

direction  of  the  mean  horizontal  velocity.  The  expression  for  Fw( k)  is  similar  to  that  for  t v(k ) 
given  the  corresponding  change  of  indices. 

Thus,  one  uses  nine  parameters  as  “input”  information  for  constructing  the  model  of  tensor  (5.2). 
The  analytical  dependences  of  /1-functions  on  these  parameters  proposed  in  (Kristensen  L.,  Lenshow 
IT  11.,  Kirkegaard  l\,  Courtney  M.,  1989)  are  not  presented  here  because  they  are  too  cumbersome, 
in  our  opinion,  an  advantage  of  (he  approach  considered  is  that  the  input  parameters  of  the  model 
can  be  obtained  in  (wo  ways:  first,  by  model  assignment  depending  on  the  conditions  of 
stratification  and  scaling  on  the  basis  of  similarity  theory  (here,  one  can  use  data  already 
published  in  order  to  reach  more  high  level  of  generalization);  second,  by  direct  determination  of 
the  unknown  parameters  on  the  basis  of  approximation  of  the  experimentally  measured  spectra  by 
expressions  (5.3)  according  to  the  least-squares  method  followed  by  calculation  of  ^-functions  and 
the  tensor  (5.2).  The  implied  possibility  of  a  model  description  of  the  spatial  structure  of 
homogeneous  anisotropic  turbulence  with  the  exclusion  of  the  necessity  to  perform  measurements  at 
many  points  is  very  attractive  irom  the  practical  point  of  view.  Hut  the  model  consideied  has  some 
shortcomings.  In  particular,  it  cannot  guarantee  positive  del initencss  of  the  tensor  (5.2)  lot 
arbitrary  input  parameters. 

The  coherence  function 


rl.  (It,  o>)  =  |  W;/R,  ©)  I V  I  Wij( 0,  to)  Wji(0,  (0)  | ,  (5.4) 

00 

whore  IV;/ 1 t,  a>)  =■•  ---  f  /f;/R,  c)exp(-i(0T)<h  is  the  mutual  spectrum  of  fluctuations  at  two 
2  K  J 

-oo 

spatially  separated  points,  is  one  of  the  most  informative  characteristics  among  others  being  used  to 
describe  statistical  laws  of  the  spatial  structure  of  the  random  wind  velocity  field. 

Let  the  observation  points  be  separated  in  a  horizontal  plane  by  the  distance  R  normally  to  the 
mean  wind  direction.  Using  the  hypothesis  of  frozen  turbulence  and  expansion  into  the  three- 
dimensional  spatial  spectrum  for  the  spatiotemporal  correlation  function  /i;/R,  t)  and  the  model 
representation  (5.2),  one  can  obtain  the  following  expressions1  for  the  mutual  spectra  of 
longitudinal,  transversal,  and  vertical  velocity  components: 

eO 

WU(R,  «)  -  2nv~]  |  /s3  k  4{k2J0(Rk)  A\(k)  +  w2  v2  [J0(Rk)  ~  R~ 1  k~ 1  J\(Rk)]A2(K)  + 
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+  to2  v  2  R  1  k  1  J \(Rk)A^(K)}dk] 


W22(R,  0))  =  2nv~'  |  k  K-'{  0)2  tr2  k2[J0(Rk)  -  R'  k~'  Jt(Rk)]A^n)  +  [(or1  vA-3R2  k2)J0(Rk)  + 

0 

+  2R- '  k(m2  V *2  +  3 R~2  )J^Rk)]A2M  +  R~2  k2[3 \J0(Rk)  +  ( R 2  k2  -  6)R~'  k~'  /,(K/e)M3(K)}d&; 

00 

Wxi(R,  w)=2nv~'  J  kK-/i{<o2v-2R-'kJt(Rk)AlM+R-2k2[3Jo(Rk)+(R2l(2-b)R-'tr'Ji(Rk)]A2(K)+ 

0 

+|<yi  -  R  2  k2)J{){RU)  -  2 (;c2  It2  -  3  R  2  k2)  R~'  k~'  Jt(Rk)\  A-jW)}dk,  (5.5) 

where  /?  -  V/;2  H-reV’2  ;  ,/„(.r)  arc  Bessel  fund  ions  of  (lie  nth  order. 

Tims  we  can  calculate  coherence,  i.e.,  obtain  information  about,  the  spatial  structure  of  an 
anisotropic  wind  velocity  field  as  a  whole  within  (lie  frames  of  the  homogeneous  turbulence 
approximation  by  using  the  parameters  of  velocity  component  autocorrelation  spectra  mcasuicd 
experimentally  at  a  single  point  as  an  input  information. 

Since  most  of  the  investigations  are  restricted  by  the  use  of  isotropic  turbulence  approximation, 
experimental  measurements  of  coherence  are  not  enough  complete  to  perform  detailed  comparison 
with  the  model  prediction. 

The  aim  of  the  study  presented  in  this  report  is  to  experimentally  verify  the  spectral  tensor  model 
of  Kristensen  L.ct  al.,  1!)8t)  which  is  not  isotropic  but.  includes  isotropy  as  a  particular  case. 

The  acoustic  weather  station,  which  is  a  compact  program-simulated  device  and  enables  one  to 
obtain  data  about  each  component  of  the  wind  velocity,  temperature,  and  their  fluctuations  at 
frequency  up  to  20  11/  and  to  measure  air  pressure  and  humidity,  was  the  principal  measuring 
instrument,  in  the  experiments  performed.  Two  to  three  such  weather  complexes  controlled  by  a 
personal  computer  were  used  in  the  experiments,  and  this  made  it  possible  to  realize  long-term 
automated  measurements  and  data  storage.  Different  spatial  arrangements  of  the  weather  stations 
both  in  the  horizontal  and  a  vertical  planes  were  used  in  different  experiments.  The  processing  of 
measurement  data  included  the  calculation  of  autocorrelation  spectra  and  coherence  spectia  for  all 
the  three  velocity  components,  friction  rate,  temperature  scale,  flows  of  momentum  and  heat, 
Obukhov  scale,  and  Richardson’s  number  what  made  it  possible  to  monitor  the  state  of  the  surface 
layer  by  stratification  conditions.  Some  preliminary  results  of  the  investigations  are  published  in 
the  1st  interim  report. 

In  general,  the  investigations  performed  confirm  the  existence  of  strong  anisotropy  of  the  wind 
velocity  field  fluctuations  in  the  surface  layer  and  the  dependence  of  the  pulsation  energy 
distribution  over  the  space  of  wave  numbers  on  the  shape  of  the  corresponding  one-dimensional 
spectra  /■„(/<■),  F„(k),  l\,Xk).  At  the  same  time,  the  /1-functions  calculated  in  accordance  with  the 
model  prove  to  be  negative  in  the  region  of  small  wave  numbers  in  some  cases  for  particular 
parameters  of  the  spectra  a2,  I,  p  obtained  experimentally.  This  leads  to  the  break  of  the  positive 
definiteness  condition  of  the  tensor  <h,;(k),  and  we  obtain  coherence  greater  than  unity  for  the 
corresponding  velocity  component. 

this  discrepancy  is  caused  by  the  drawbacks  of  the  model  related  to  the  fact  that,  in  real 
experiments,  one  often  observes  nonslationarity  caused  by  daily  radiation  fluctuations  and  leading 
to  change  of  the  heat  flow  II  depending  on  time  and  height.  This  effect  is  especially  strong  under 
the  spacing  R  considerably  exceeding  the  value  of  the  outer  turbulence  scale. 
I  low  ever,  the  contribution  of  negative  values  ol  /13  corresponding  to  liirgc-sizc  vortices  decreases  fot 
horizontal  spacing  equal  or  somewhat  less  than  /„,  and  the  model  coherence  does  not  exceed  unity. 
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The  spectra  of  wind  velocity  components  measured  over  the  plain  underlying  surface  at  midday  time 
in  summer  under  unstable  stratification  of  the  surface  layer  for  z/L  =  -0  064  and  normalized  by 
(Ik-  variance  are  presented,  as  an  example;  in  Fig.  13.  The  sample  contains  8192  synchronously 
measured  values  of  three  wind  velocity  components  and  temperature.  the  measurements  wue 
performed  at  a  rate  of  10  11/  during  13.65  min. 


/.  /<•;.  / One  dimensional  spectra  for  the  longitudinal,  transverse,  and  vertical  wind  velocity 
components ,  u ,  v,  w,  obtained  experimentally. 


riG.  PI.  Three  model  A  functions  determining  the  spectral  tensor,  calculated  for  the  spectra 

depicted  in  Fig.  13- 

The  whole  sample  was  divided  into  16  sequential  intervals  containing  512  potato  each  ir >  onferto 
increase  the  degree  of  freedom  in  the  spectral  processing.  Linear  filtration  of  data  was  earned  out. 
using  the  approximation  with  a  polynomial  of  the  first  degree  followed  by  subtraction  from  m 
initial  data  array  So  the  trend  effects  in  the  wind  velocity  were  excluded.  1  he  spectra  calculated 
for  each  interval  were  summed  up  and  normalized  by  the  sample  variance  Hie  resuJ ted  spectra 
were  smoothed  over  five  neighboring  points.  According  to  our  estnnat^ 
spectrum  estimation  using  this  processing  does  not  exceed  1  A. 
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components  averaged  over  measurement  time  were  as  follows:  Vu-  1.446  m/s,  V„  -  0.000  m/s, 
y  =  -0.086  m/s.  Spectral  parameters  obtained  by  the  least-squares  method  for  these  velocities 
are  equal  to 

al  =  0.190  m2/s2  al  =  0.145  m2/s2  al  =  0.056  m2/s2 

/„  =  1.074  m  lv  =  1  095  in  lw  =  0.426  m 

p„  =1.2  Po  ~  1  -5 

Figure  14  presents  three  /1-functions  determining  the  spectral  tensor  (5.2)  and  calculated  for  the 
above-mentioned  parameters  in  correspondence  with  the  model  of  Kristcnsen  L.et  al.,  1989. 
Figure  15  presents  the  coherence  spectra  calculated  by  formulas  (5.4)  and  (.>.5),  and  those  obtained 
experimentally  at.  the  transverse  horizontal  separation  ft  —  0.5  m.  A  little  excess  of  the 
experimental  coherence  over  the  calculated  one  can  lie  explained,  on  the  one  hand,  by  the 
shot  I  comings  of  the  model  restricted  by  the  hypotheses  of  frozen  and  homogeneous  turbulence  and, 
oa  (he  olli'T  hand,  by  errors  in  determining  input  parameters  of  the  model  and  in  performing  the 
mode]  calculations.  Since  the  model  considered  imposes  increased  requirements  on  the  correctness 
of  slat  is*  Seal  processing  of  the  measurement  data,  it  is  necessary  to  eliminate  low-frequency  trends 
(In'  initial  temporal  series  of  data  in  order  to  increase  the  reliability  of  spectra!  estimations  in  the 
of  low  frwjiH  iicp'S. 


I'lC.  I1).  Coherence  spectra  for  longitudinal  (I,  1'),  transverse  (2,  2’),  and  vertical  (3,  3') 
wind  velocity  components  obtained  experimentally  (curves  1~3)  and  calculated  by  the  mod  cl  (!'- 


V). 


On  (he  basis  of  a  good  coincidence  of  the  predicted  and  the  experimentally  measured  coherence,  one 
can  come  to  a  conclusion  that  the  model  considered  can  be  useful  in  estimating  spatial  statistics  of 
turbulent  pulsations  using  information  obtained  at  a  single  point  of  the  space  when  the  conditions 
of  the  experiment  correspond  to  the  accepted  hypotheses. 
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6  STUDY  OF  STATISTICAL  PROPERTIES  OF  THE  WIND  AND 
TEMPERATURE  FIELDS  IN  THE  ATMOSPHERIC  SURFACE  LAYER 
WITH  ULTRASONIC  SENSORS 


At  present  to  develop  and  substantiate  theoretically  the  lidar  and  acoustic  methods  ot  sounding 
the  atmospheric  parameters  normally  the  Kolmogorov  model  of  turbulence  for  locally  isotropic 
medium  is  used.  In  this  ease  a  spatiotemporal  structure  of  wind  velocity  field  is  described  using 
the  hypothesis  of  “frozen”  turbulence,  and  temperature  is  considered  as  a  conservative  and 
passive  impurity.  But,  as  known,  in  the  atmospheric  surface  layer  the  turbulence  is  essentially 
anisotropic.  Moreover,  there  are  differences  in  the  behavior  of  the  wind  and  temperature 
fluctuations  depending  on  the  stability  state  of  the  atmosphere  (Newstad  F.T.M.  and  Van  Dop 
II.,  1985) 

In  this  report  we  consider  some  results  of  a  series  of  experiments  to  study  in  detail  the  structure 
of  turbulence  in  the  atmospheric  surface  layer  using  six  ultrasonic  sensors  of  wind  velocity  and 
temperature  fluctuations.  The  primary  goal  of  the  study  was  to  elucidate  a  dependence  of 
statistical  characteristics  of  fluctuations  of  the  wind  velocity  and  temperature  on  a  variety  ot 
parameters:  the  distance  between  the  measurement  points,  direction  of  separation  with  respect 
to  average  velocity,  frequency,  and  thermal  stratification  conditions  in  the  surface  layer. 


Measurements  were  carried  out  during  twenty  four  hours  every  half  an  hour.  Spatially  the 
ultrasonic  sensors  were  placed  along  one  line  at  a  distance  of  1  m  between  them.  Measurements 
of  fluctuations  of  the  longitudinal  wind  velocity  component  and  temperature  were  made  with 
all  the  six  sensors.  Moreover,  the  acoustic  weather  station  (Rostov  A. I'.,  1993  ),  which  measured 
(lie  fluctuations  of  three  wind  velocity  components,  temperature,  humidity,  and  average  value 
of  air  pressure,  was  used  as  one  of  the  sensors.  The  measurements  were  made  synchronously  at 
a  rate  of  4  llz.  The  length  of  a  single  data  sample  was  4096  points,  i.e.  the  averaging  time 
equaled  about  17.1  min. 

The  study  was  carried  out  in  three  successive  stages  in  accordance  with  the  longitudinal, 
vertical,  and  transverse  disposition  of  the  sensor  linearly  arranged  set,  relative  to  the  average 
wind  direction.  At  each  stage  the  experiment  was  carried  out  in  an  automated  mode  o 
operation  controlled  using  two  personal  computers.  These  computers  also  recorded  and  partially 
processed  the  data  on  a  real  time  scale.  Thus  recorded  data  have  then  been  subjected  to  a  more 
detailed  analysis  and  statistical  processing. 

We  have  analyzed  35  measurement,  sessions.  For  each  run  the  following  statistical 
characteristics  were  determined:  for  the  longitudinal  wind  velocity  component  and  temperature 
we  calculated  the  average  values,  variances,  asymmetry  and  excess,  autocorrelation  functions, 
and  auto  spectra  at  each  sensor  point,  and  also  the  mutual  correlation  functions,  spectra  of 
coherence  and  phase  were  calculated  for  the  sensors  placed  at  the  distance  1  5  m  from  each 
other  Besides,  based  on  the  data  of  the  weather  station  we  calculated  the  single-point  moments 
and  spectra  for  the  transverse  and  vertical  wind  velocity  components  and  the  scales  o  velocity 
(friction  velocity),  temperature,  specific  humidity,  length  (Obukhov  scale),  vertical  turbulent 
flux  of  the  momentum,  water  vapor,  explicit  and  latent  heat  were  also  calculated. 


Although  the  question  on  similarity  of  temperature  fluctuations  to  the  wind  ones  has  been 
widely  discussed  by  Byzova  N.L.  ,  Ivanov  V.N.  ,  and  Garger  E.K.,  198.))  at  present  an 
assumption  about  the  identical  nature  of  the  turbulent  heat  and  momentum  exchange  in  ie 
surface  layer  of  the  atmosphere  is  a  hypothesis  and  needs  for  the  further  experimental  tests.  In 
this  paper  we,  however,  are  limited  by  a  comparative  analysis  of  the  behavior  of  spatiotemporal 
correlation  functions  and  coherence  spectra  of  fluctuations  of  the  temperature  and  wind  velocity 
components. 


The  range  of  variation  of  the  thermal  stratification  conditions  during  the  experiment,  in  terms 
of  the  parameter  z/L  (where  /.  is  the  Obukhov  scale,  2  is  the  observation  altitude),  was  quite 
broad,  i.e.,  the  stale  of  the  surface  layer  changed  from  a  very  stable  state  (1  class,  by  lemu) 


at  night  to  a  very  unstable  one  (7  class)  during  day-time.  Change  of  the  surface  layer  state  is 
shown  in  Fig.  16. 

As  known  (Lumley  J.L.  and  Panofsky  II. A.,  1964  ),  the  temperature  fluctuations  arise  in  the 
atmosphere  only  if  the  height  stratification  of  the  potential  temperature  occurs.  According  to 
(Lumley  J.L.  and  Panofsky  II. A.,  1964  )  no  temperature  fluctuations  are  generated  in  the 
presence  of  adiabatic  gradient  of  the  potential  temperature  (neutral  stratification  conditions). 
When  the  vertical  temperature  gradient  deviates  from  the  adiabatic  one,  vertical  motion  of  the 
elementary  volumes  of  air  by  fluctuations  of  the  vertical  velocity  leads  to  local  temperature 
fluctuations.  Reduction  of  the  variance  of  these  fluctuations  a]  to  a  dimensionless  form  is 

done  using  the  temperature  scale  T*  determined  by  the  equality  <t'w'>  =  ~u*T*,  where  V  and 
w*  are  (he  fluctuations  of  temperature  and  wind  velocity  vertical  component,  u*  is  the  friction 
velocity. 


FIG.  16. 


Figure  17  shows  the  normalized  variance  of  temperature  fluctuations  as  a  function  of  stability 
according  to  data  from  our  measurements.  One  can  see  that  at  the  stratification  close  to  neutral 
the  magnitude  at/  \  T*  |  increases.  On  the  whole  the  dependence  presented  well  agrees  with  the 
theoretical  statements  and  known  experimental  data. 
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FIG.  17. 


In  (Byzova  N.L.  ,  Ivanov  V.N.  ,  and  Garger  H.K.,  1989  )  the  similarity  principle  of  temperature 
fluctuations  spectra  to  wind  ones  is  used.  However,  there  is  a  difference  between  these  spectra 
connected  with  the  temperature  field  may  not  be  stretched  itself,  but  it  is  transferred  by  the 
velocity  field.  In  this  case  the  inertia  intervals  of  these  two  fields  overlap  partially,  but  their 
coincidence?  is  not  necessary  (Lumley  J.L.  and  Panofsky  H.A.,  1964  ).  A  possible  cause  of  the 
difference  in  temperature  and  wind  velocity  spectra  may  be  different  behaviors  of  the  turbulent 
exchange  coefficients.  At  present  the  difference  between  the  turbulent  exchange  coefficients 
depending  on  the  stratification  is  established  quite  reliably  (Newstad  F.T.M.  and  Van  Dop  IL, 
J9"5).  In  accordance  with  the  Mon  in -Obukhov  similarity  theory  in  the  atmospheric  surface 
layer  the  ratio  between  the  turbulent  exchange  coefficients  a  for  the  heat  /Cj,  and  momentum 
/sTin  is  a  function  of  the  stability  parameter  (dimensionless  height) 
£  =■ z/L  :  a  =  K\}/Km  =:  <Pm(£)/<Ph(£),  but.  it  is  impossible  to  determine  the  functions  (pm(£) 
and  <p| ,(^)  with  the  similarity  theory.  To  do  this  one  needs  for  corresponding  experimental  data. 
In  accordance  with  the  results  of  numerous  measurements  of  the  velocity  and  temperature 
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profiles,  including  ( he:  laboratory  conditions,  the  value  a  was  found  to  be  1.16  at  £  =  0.  Since, 
by  definition,  cpul(0)  -  1  the  function  q>|,(0)  -  0.86.  However,  many  authors  consider  that  the 
relations  of  q>m  and  to  £,  obtained  by  processing  of  the  experimental  data  requires  further 
refinement  (Newslad  F.T.M.  and  Van  Dop  II.,  1985). 

Mxaiuples  of  the  temperature  and  wind  velocity  fluctuation  spectra  normalised  to  the  variance 
we  have  obtained  under  conditions  of  stable,  neutral,  and  unstable  stratification  are  presented 
in  Figs.  18,  19,  and  20,  respectively.  Unfortunately  for  the  stable  and  neutral  states  of  the 
surface  layer  stratification  the  intensity  of  turbulent  fluctuations  was  very  weak,  and  the 
spectral  characteristics  were  close  to  those  of  pure  noise.  At  large  z/L  values  the  short- 
duration  sporadic  manifestations  of  the  turbulence  lasting  0.5-5  min  were  observed  in  certain 
cases  only.  Since  the  averaging  time  was  about  17  min  this  led  to  a  disturbance  of  the 
stationary  conditions  during  the  sampling  time.  It  is  clear  from  Figs.  18  and  19  that:  the  inertia 
interval  in  spectra  is  small,  and  a  rise  in  the  spectral  curves  in  the  high  frequency  region  is 
explained  by  noise.  These  results  agree  with  the  conclusions  from  Byzova  N.L.  ,  Ivanov  V.N.  , 
and  Garger  F.K.,  1989  that  in  the  night  surface  layer  the  turbulence  is  intermittent.  Also  these 
results  agree  with  I  he  data  of  similar  experiments  described  by  Newslad  F.T.M.  and  Van  Dop 
IF,  1985,  where  il  was  shown  that  at  Ri  ~  0.2  in  some  observational  series  the  slope  ol  the 
dimensionless  spectra  multiplied  by  the  frequency  was  close  t;o  zero  and  had  the  values  from  - 
1.2  to  2.5  in  the  other  ones. 


Measurements  of  temperature  fluctuations  under  unstable  stratification  showed  that  almost  all 
disturbance  energy  concentrates  near  the  spectral  maximum  that  is  clearly  observed  in  the  low 
frequency  region,  and  in  the  inertia  interval  the  spectra  really  obey  the  “2  thirds  law”  (see 
Fig.  20). 
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n  is  known  (Ncwslad  l\T.M.  and  Van  Hop  II.,  1085,  Byzova  N.L.  ,  Ivanov  V.N.  ,  and  Gnrgcr 
i*.K.,  1989,  I  aim  ley  J.l  .  and  Panofsky  1 1. A.,  1964  )  dial  in  the  surface  layer  the  turbulence  is 
essentially  anisotropic.  I, el.  us  consider  some  manifestations  of  the  anisotropy  in  the  example  of 
spat  iol  emporal  correlal  ion  functions  of  find  uat  ions  of  Hu*  longitudinal  wind  velocity  component 
and  temperature?  measured  under  the  condition  of  unstable  stratification,  which  are  shown  in 
Fii's-  21 ,  22,  and  2d  for  the  longitudinal,  transverse,  and  vertical  separations  between  the  sensor 
points  of  the  acoustic  meters,  respectively. 

The  most  slow  drop  of  the  correlation  maxima  takes  place  in  the  longitudinal  direction  both  for 
the  velocity  fluctuations  (the  curves  /-•/)  and  for  temperature  ones  (the  curves  /W')- 
Moreover  when  the  spatial  separation  of  points  increases  (the  curves  /  to  4  correspond  to  the 
separations  1,  3,  4,  and  3  m),  the  natural  decrease  of  the  correlation  maxima  and  rise  in  the 
time  delay  of  their  location  are  observed,  i.e.  a  satisfactory  fulfillment  of  the  hypothesis  of 
“frozen”  turbulence  takes  place  in  this  direction. 
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In  the  transverse  direction  the  drop  of  the  wind  fluctuation  correlation  occurs  considerably  more 
quickly  then  in  the  longitudinal  direction,  and  the  level  of  the  temperature  correlation  is 
considerably  lower  then  for  wind.  Moreover,  if  for  fluctuations  of  the  longitudinal  velocity  a 
correspondence  of  the  maximum  values  and  their  temporal  delays  to  the  values  of  spatial 
separations  is  still  traced,  then  for  the  temperature  fluctuations  this  dependence  is  broken,  and 
a  shift  of  maxima  is  absent  at  all. 

As  is  clear  from  Fig.  23  for  the  separation  of  meters  in  the  vertical  direction  (here  the  curves  1- 
4  correspond  to  the  separations  1,  2,  3,  and  4  m)  the  correlation  of  the  longitudinal  velocity 
fluctuations  was  most  low.  However,  here  a  decrease  of  the  statistical  relation  of  fluctuations  is 
traced  as  far  as  the  difference  of  heights  increases.  The  peculiarity  of  the  fluctuation  correlation 
with  the  vertical  separation  is  that  the  correlation  maximum  between  the  levels  (for  example, 
the  levels  at  2  and  3  m  apart  -  the  curves  2  and  3,  respectively)  is  observed  not  for  the  positive 
time  delay,  as  for  the  cases  of  the  longitudinal  and  transverse  separations,  but  for  the  negative 
ones,  i.e.  when  the  observations  at  the  upper  level  precede  to  observations  at  the  lower  level. 
Such  a  shift  of  the  cross-correlation  function  maximum  in  accordance  to  Lumley  J.L.  and 
Panofsky  H.A.,  1964  means  a  slope  of  turbulent  eddies  towards  the  average  motion  due  to  their 
interaction  with  the  positive  (upward)  gradient  or  the  average  velocity. 

It  is  interesting,  that  when  the  distance  between  the  levels  is  not  too  large  (R  =  1  m,  the  curve 
/)  in  the  given  realization  a  presence  of  two  weakly  pronounced  local  maxima  with  the  time 
delay  opposite  in  sign  is  observed.  It  is  possible,  that  the  right-hand  side  maximum  with  the 
positive  delay  is  connected  with  the  vertical  component  of  the  average  wind  velocity  caused  by 
the  air  flows  from  heated  underlying  surface  occurs  in  the  surface  layer  under  conditions  of 
unstable  stratification.  It  must  be  so,  that  the  regular  vertical  transfer  of  inhomogeneities  by 
this  component  leads  to  appearance  of  maximum  at  a  positive  delay.  Since  this  component  is 
small  (in  this  case  <w>  =  0.1  m/s)  its  influence  on  the  maximum  shift  is  detected  only  for 
close  levels,  where  the  influence  of  corresponding  difference  of  the  longitudinal  velocities  caused 
by  the  presence  of  the  logarithmic  profile  is  less  important .  As  one  can  see  from  Fig.  2o  the 
right-hand  side  maximum  practically  disappears  at  the  separation  R  =  2  m  already  (owing  to 


27 


the  influence  of  ordinary  processes  of  the  evolution  of  inhomogeneities  of  the  velocity  and 
turbulent  mixing). 

As  regards  the  temperature  fluctuations  although  the  correlation  is  also  low,  as  in  the  case  of 
the  transverse  separation,  one  can  notice,  that  the  maximum  of  functions  is  reached  not  at 
negative  time  shift,  as  for  the  longitudinal  velocity  fluctuations,  but  for  a  positive  one  i.e.  y 
the  transfer  of  temperature  inhomogeneities  upwards  by  the  average  vertical  component  of  the 

wind  velocity. 

On  the  whole  Figures  21-23  show  that  during  the  measurements  under  moderately  unstable 
stratification  conditions  of  the  surface  layer  the  velocity  vortexes  were  most  strongly  stretched 
along  the  direction  of  average  motion.  Typical  dimension  of  the  vortexes  in  the  transverse 
direction  was  smaller  than  in  the  longitudinal  one,  but  it  was  considerably  larger  than  in  the 
vertical  direction.  Spatial  correlation  of  the  temperature  inhomogeneities  drops  most  slowly  tor 
the  horizontal  position  of  the  sensors  linear  array  along  the  average  wind  direction  also_ 
However  in  the  transverse  and  vertical  directions  the  behavior  of  the  correlation  functions 
the  temperature  fluctuations  essentially  differs  from  that  of  corresponding  corre  ation  functions 
of  the  longitudinal  velocity  fluctuations.  All  this  shows  that  the  essential  distinctions  in  the 
structure  of  turbulent  fields  of  wind  velocity  and  temperature  which  are  expressed  in  different 
manifestations  of  the  spatial  anisotropy  of  their  statistical  characteristics. 


To  study  the  detailed  structure  of  turbulent  fields  of  the  velocity  and  temperature  information 
on  the  degree  of  statistical  relation  of  fluctuations  in  different  spatial  points  for  the  specific 
values  of  frequency  is  of  certain  interest  .  Such  informat  ion  is  m  the  coherence  functions 


Coh'V.  «;>  =  l  Co2(/)  +  Q2(P  1  /  I  f'aW  F‘b(PI  i-  )  ~  C  2-  3- 


(6.1) 


where  f  is  the  frequency;  F\(P  and  F‘B(P  are  the  power  spectra  at  the  points  A  and  B, 

respectively  R;  is  the  distance  between  the  points  A  and  B;  Co (/)  is  the  cospectrum,  and  Q(p 
is  the  quadrature  spectrum;  indices  i,  ;  =  1.  2,  3  mean  the  longitudinal,  transversal,  and  vertical 
wind  components  and  longitudinal,  transversal,  and  vertical  disposition  of  the  sensor  array 
relative  to  the  average  wind  direction,  respectively. 

The  coherence  (6.1)  is  a  statistical  characteristic  determining  the  contribution  of  fluctuations  of 
different  frequencies  into  the  correlation  between  the  fluctuations  of  meteorological  quantities 
at  two  spatially  separated  points.  As  applied  to  the  wind  velocity  the  coherence  can  be  treated 
as  a  measure  of  the  velocity  vortex  stability. 

Figures  24  and  25  show  the  behavior  of  the  coherence  functions  for  the  horizontal  velocity 
components  (curves  1  and  2)  and  for  temperature  (curves  V  and  20  measured  for  the 
separations  between  the  observation  points  R  =  1  m  (curves  /  and  1')  and  R  5  m  (curves  2 
and  20  along  the  average  wind  velocity  direction  for  conditions  of  unstable  and  stable 
stratifications,  respectively.  One  can  see  from  Fig.  24  high  values  of  the  coherence  between 
fluctuations  at  the  points  of  sensor  locations  in  the  low  frequency  region  for  unstable 
stratification,  as  it  could  be  expected  from  the  known  conceptions  of  the  turbulent  motion  in 
the  near  surface  layer  (Lumley  J.L.  and  Panofsky  H. A. ,1964). 


When  the  frequency  increases,  i.e.  when  going  to  smaller  scales  of  turbulent  inhomogeneities,  a 
decrease  in  the  coherence  takes  place  as  in  the  case  when  a  distance  between  the  measurement 
points  increases.  In  this  case  obvious  distinctions  in  the  behavior  of  temperature  and  wind 
fluctuations  are  not  revealed  in  the  case  when  the  correlation  functions  at  the  longitudinal 
separation  of  sensors. 

When  the  stratification  is  stable  (Fig.  25),  a  drop  of  the  coherence  with  increasing  frequency 
and  distance  between  the  points  is  more  quick  than  it  occurs  under  unstable  one.  This  tact  can 
be  explained  by  a  general  decrease  in  the  turbulence  energy  due  to  the  buoyancy  Archemedes 
force  and  corresponding  decrease  of  all  scales  (Newstad  F.T.M.  and  Van  Dop  H.,  1985). 
Moreover,  a  certain  distinction  in  the  behavior  of  temperature  and  wind  coherence  is  observed 
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the  typical  scale  of  the  longitudinal  velocity. 

If  one  takes  into  account  that  the  mechanism  of  temperature  fluctuation  generation  for  the 
stratification  different  from  neutral  is,  on  the  whole,  caused  by  the  vertical  velocity  component 
fluctuations  then  such  a  scale  must  be  certain  typical  size  of  the  turbuient  inhomogenei  les 
along  the  vertical  direction.  Under  stable  stratification  conditions  the  turbulent  eddies 
velocity  are  compressed  along  the  vertical  direction  by  the  buoyancy  force  , 

Thus  the  aforesaid  distinctions  of  the  turbulent  exchange  coefficients  and  spatial  anisotropy  of 
th^turbulen^hnhomogeneities  is  also  a  factor  causing  the  difference  in  the  behavior  of  spectral 
characteristics  of  the  longitudinal  velocity  component  and  temperature  in  the  low  frequ  y 
Son ^  At  higher  frequency  time  scales  are  smaller  and  the  eddies  are  more  isotropic. 

For  a  theoretical  generalization  of  the  experimental  data  of  wind  coherence  the  following  model 
representations  (Mazurin  V.P.,  1983)  are  often  use,4 

Coh’V,  Rj)  =  exp  (-  K 'j  rij)  i,  j  =  1,2,3,  (6-2) 

where  nj-fRj/Vi  is  the  dimensionless  frequency;  K)  is  the  attenuation  parameter.  The 
,  %  in  Fine  24  and  25  are  the  exponential  approximation  of  the  corresponding 

performed  by  different  authors  in  different  sites. 

According  to  the  conclusions  from  Mazurin  V.P.,  1983,  under  unstable  stratification  no 
noticeable  dependence  of  k\  on  R  exists.  Really,  the  values  of  attenuation  parameter  Kx  f 
i?=1  and  5  m  obtained  from  our  measurements  have  close  magnitudes  of  5.03  and  4  01 
respectively  (the  curves  la  and  in  Fig.  24).  At  the  same  time,  under  stable  stratification 
(curves  la  and  2a  in  Fig.  25)  these  parameters  are  4.92  and  10.1  for  the  same  separations. 
These  values  well  agree  with  the  conclusions  from  Mazurin  V.P.,  1983  that  the  growth  of  Kx  is 
observed  at  the  increase  of  R  under  stable  stratification.  Moreover,  if  at  ^  »eparations 
(R  =  5)  the  growth  of  k\  is  observed  when  the  stability  increases  (the  values  4.92  and  10.1  tor 

the  unstable  and  stable  conditions),  at  small  separations  (R  =  1)  the  attenuation  parameter  can 
even  decrease  (the  values  5.03  and  4.01,  respectively).  Such  a  dependence  agrees  with  the 
results  of  some  investigations  described  in  Mazurin  V.P.,  1383. 
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FIG.  25 

To  obtain  a  more  detailed  quantitative  conclusions  about  the  influence  of  the  stability 
conditions  in  the  surface  layer  on  the  degree  of  the  turbulence  spatial  anisotropy  more  long 
experiments  are  needed  allowing  for  the  influence  of  a  set  of  additional  factors  like  the 
boundary  layer  height,  underlying  surface  roughness,  etc.  Nevertheless,  it  is  obvious  from  the 
analysis  presented  that  similar  complex  automated  measurements  are  a  good  possibility  to  study 
the  surface  layer  turbulence  where  the  synchronous  data  sampling  is  needed  for  a  large  number 
of  hydrometeorological  parameters.  The  cases  of  weak  and  intermittent  turbulence*  require 
additional  study.  A  positive  aspect  applying  the  acoustic  meters'*  from  this  point  of  view  is  high 
accuracy  and  high  rate  of  recording  velocity  and  temperature  fluctuations  allowing  the  tine 
effects  in  structure  of  turbulent  fields  to  be  investigated.  The  data  obtained  can  be  used  for 
theoretical  generalization  and  identification  of  new  relationships  of  the  behavior  of  small-scale 
turbulence  in  the  atmospheric  surface  layer. 
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7. STRUCTURE  CHARACTERISTIC  FOR  TEMPERATURE  FLUCTUATIONS 
AND  THE  OUTER  SCALE  OF  ATMOSPHERIC  TURBULENCE 
RECONSTRUCTED  FROM  THE  DATA  OF  ACOUSTIC  SOUNDING 


Nowadays  acoustic  radars  (sodars)  are  widely  used  to  measure  the  profiles  of  structure 
characteristics  for  temperature  fluctuations  in  the  atmospheric  boundary  layer.  However  a 
comparison  of  sodar  data  with  in  situ  measurements  shows  that  the  standard  deviation  of  sodar  data 
reaches  40%  at  altitudes  30-500  m.  Their  maximum  discrepancy  may  reach  a  factor  of  2-5.  We 
believe  that  a  possible  reason  for  the  discrepancy  is  neglect  of  the  turbulent  attenuation  of  a  sodar 
signal  when  it  propagates  to  a  sounded  volume  and  back  to  a  receiver.  In  the  present  report  the 
interpretation  of  sodar  data  is  performed  with  the  use  of  an  iterative  algorithm  (Krasnenko  N.P 
Shamanaeva  L.G.,  1997)  considering  the  turbulent  signal  attenuation.  This  allows  the  vertical 
profiles  of  the  outer  scale  of  atmospheric  turbulence  to  be  obtained. 

The  outer  scale  of  turbulence  plays  an  important  role  in  the  theory  of  turbulence,  because  it  defines 
the  low-frequency  boundary  of  the  inertial  subrange  in  the  spectra  of  temperature  and  wind  velocity 
fluctuations  (Kolmogorov  A.N.,  1941).  Knowledge  of  L0  is  necessary  for  estimating  the  turbulent 
attenuation  of  a  sound  wave  (Baikalova  R.A.  et  al  1988).  Nevertheless  the  vertical  profiles  of  JLo 
have  never  been  measured  in  the  atmospheric  boundary  layer,  although  Sazarm  M.,  1990  poin 
their  strong  dependence  on  local  orography. 

In  the  present  report,  we  discuss  the  vertical  profiles  of  the  structure  characteristic  for  temperature 
fluctuations  c\  and  the  outer  scale  of  atmospheric  turbulence  L0  measured  with  the  Zvuk-2  three- 
channel  Doppler  sodar.  Sodar  specifications  are  given  below. 


Vertical  sensing  range,  m 
Operating  frequency  f,  Hz 


....  40-700 
1650-1850 


Electrical  transmitted  power  P0,  W . 

Pulse  repetition  period  At,  s  . . 

Pulse  duration  x,  ms  . 

Effective  transceiving  antenna  area  Ac ff,  m2 
Efficiency  of  conversion  of  electrical  power 
into  acoustic  power  and  vice  versa,  yi  =  72 


.  65 

4,6,  12 

.  150 

.  0.5 

.  0.1 


The  structure  characteristic  of  temperature  fluctuations  is  calculated  by  the  formula 


,,;w  v  2.7  •  102  P(zi)z]  (7.1) 

CfJ'  (*j)  =  - rO-bTT’ 

Y,  y2  PoCAefrrt  ( Zi ) 

where  j  is  the  iteration  number,  j  =  1 . M;  z;  is  the  current  altitude  of  the  sounded  volume, 

2-  =  iAz,  Az  =  cx/2  is  the  altitude  resolution  of  the  sodar,  i  -  i0,  to  +  1,—.  *0  specifies  the 

minimum  sounding  altitude  z0  =  i0Az,  N  specifies  the  maximum  sounding  altitude  and  is  equal  to 

the  number  of  range  gates  in  the  realization;  c  is  the  sound  speed  in  air,  c  =  20.05  Vr  ;  T  is  the 
surface  air  temperature;  X  =  c/f  is  the  acoustic  radiation  wavelength;  P(z,)  is  the  power  of  the 
acoustic  signal  coming  from  the  altitude  Z{\ 

=  L.M  Ltulizi)  (7’2) 

is  the  attenuation  factor  considering  the  signal  attenuation  on  the  propagation  path  to  the  sounded 
volume  and  back  to  the  receiver  caused  by  the  classical  and  molecular  absorption  Labs(2i>  and  by 

the  turbulent  attenuation  L  nil  (z*)  due  to  scattering  by  temperature  and  wind  velocity  fluctuations. 

The  first  iteration  was  calculated  from  Eq.  (7.1)  without  turbulent  attenuation  for 
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L(0)(z,)  =  Labs<z,)  =  exp  {-  2  (pc,  +  pmo|)  *}•  (7  ^ 

In  so  doing,  the  classical  and  molecular  absorption  coefficients  were  calculated  by  the  formulas 
(Neff  W.D.,’  1975) 


pc, (/)  =  4.0210  11  f2,  [m  ']  , 


(7.4) 


PmolV >  t,  u)  = 


Pa 


304.8 


(0.1 8^)2 


v\  ^ 


1  +  /a 


1/2 


,1m']. 


(7.5) 


where  fy  =  f/fm,  fm  is  the  frequency  at  which  the  maximum  absorption  occurs  at  the  given 
temperature  t[°C],  pressure  p[mbar],  and  relative  humidity  m(%],  determined  from  the  formula 

fm  =  (10  +  6600  h  +  444000  h 2)  p\/ ■  (76) 

Here,  the  normalized  pressure  P\=p/ 1014,  tt  =  (1.8/  +  492)/519,  and  the  parameter 
h  =  evm/px  depends  on  the  water  vapor  pressure  e,  in  the  atmosphere  which  is  given  by  the 
formula  (Komarov  V.S.,  1971) 


e,  =  0.061  1107  10M,m  (77) 

with 

Mi  =  7.63  t/(241.9  +  t)  . 

Then  the  parameter  pmax  in  Eq.  (7.5)  is  calculated  from  the  relation 

pmax  =  0.0018  tf  5  exp  {7.77  (1  -  tj'))  -  (79) 

The  second  and  the  next  iterations  of  c\  were  calculated  considering  the  turbulent  attenuation. 
The  turbulent  attenuation  factor  in  Eq.  (7.2)  (for  the  ;th  iteration)  was  calculated  by  the  formula 


L  turl  (z,)  =  exP|-2Az|;k>(2/)  +  ^^)]}  ■  (710) 

Here,  the  coefficients  of  scattering  by  the  temperature  and  wind  velocity  fluctuations  were 
calculated  by  the  formulas 


flfy'Hzi)  =  9.545  c7v[Zi) 

c  X 

ftp  (z,)  =  0.3596  fa— 

r'  T  X 


(7.11) 

(7.12) 


derived  in  (Baikalova  R.A.  et  al„  1988)  for  l\  » which  is  satisfied  in  our  case  (X  -  20  cm). 
The  ;th  iteration  of  c\  entering  Eq.  (7.12)  was  calculated  by  Eq.  (7.1). 

The  outer  scale  of  atmospheric  turbulence  and  the  structure  characteristic  of  wind  velocity 
fluctuations  were  determined  as  follows. 

The  structure  characteristics  of  the  wind  velocity  fluctuations  were  determined  by  averaging  of  the 
instantaneous  wind  velocity  components  VXfc(zf),  V^(z,),  and  V^(z,)  measured  by  the  Zvuk-2 

sodar 
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K- 1 


cl  (z;)  =  (Atjviyiz/) + v2y*Az,)+  vi,Azi))  /  £  [^i)  -  vk+iUi)]2 


(7.13) 


k  =  1 


where  =  12  s  is  the  sound  pulse  repetition  period,  M  -  50  is  the  number  of  pulses  in  a 
measurement  series, 


K  !  A 

vxM=-^  £  vXk(Zi),vyM=T  £  vyk(zi), 

A  ,  .  r  _  j 


k  =  1 


k  =  1 


K 


£  v,kUi) 


k  =  i 


are 


the  components  of  the  wind  velocity  averaged  over  the  period  KM,  and 
Vxk  (2/)  ■  Vx  av  ( z, )  +  Vyk  (zj)  ‘ Vyav  izj)  +  Vjk  izj)  '  V  ?  w  (fj) 


Vk(zi)  = 


ajv  iZf)  ^  V yav  (zi)  +  V zav  ( 2/ ) 


(7.14) 


(7.15) 


Because  in  our  experiment  we  measured  only  the  wind  velocity  components  averaged  over  10-min 
intervals,  d  was  approximated  by  the  relation 

Cl  (*i)  =  CVo  (0.04  +  0.33  zp  )  (716) 

suggested  in  (Brown  E.H.and  Clifford  S.F.,  1976)  for  the  daytime  convection.  The  parameter  CV(> 
was  determined  from  the  data  of  in  situ  measurements  performed  by  micropulsation  sensors  placed 
17  m  above  the  ground. 

The  outer  scale  of  atmospheric  turbulence  was  estimated  as  follows.  At  visible  wavelengths, 
Tatarskii  V.I.,  1961  derived  the  expression  for  the  structure  characteristic  of  the  refractive  index 


d (^f)  =  2.8  M(z)  L4JHz), 


(7.17) 


where  M(z)  is  the  potential  refractive  index  gradient  related  with  the  gradient  of  the  potential 
temperature  logarithm  by  the  formula 


Miz)  =  -  78  •  10 


,6  £(f)  Sin 6(a) 
T(z )  52 


(7.18) 


In  its  turn,  the  potential  temperature  is  defined  by  the  air  temperature  and  pressure 


0(2)  =  T(z)  [1000/p(2)]°'286. 


(7.19) 


According  to  the  formula  (Coulman  C.E.  et  al.,  1988) 

c5(2)  =  CfW'78-10^2 


Tl(z) 


(7.20) 


the  structure  characteristic  of  the  refractive  index  at  visible  wavelengths  is  connected  with  the 
structure  characteristic  of  the  temperature  fluctuations,  so  from  Eqs.  (7. 17)— (7. 19)  we  derive 


L0(z)  = 


C\(z) 

2(50/ 5  2]2 

1,1000  J 

3/4 


(7.21) 
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For  the  atmospheric  boundary  layer  at  altitudes  40  S  *  <  700  m,  we  can  take  p(z)  =  *00° 

50/82  =  ya  (see  Matveev  L.T.,  1976),  where  ya  is  the  adiabatic  temperature  gradient.  Fi  y, 
the  ;th  iteration  of  the  outer  scale  of  atmospheric  turbulence  we  obtain 

(7.22) 


4>)(^)  =  [C^)Oi)/(2.8^)]3/4. 

The  iterative  algorithm  of  acoustic  sounding  data  processing  was 

The  structure  character, ich’sSc  ^  f“o. 
scale  o  a^”os^  er“.ni"lim  sounding  altitude  r„  -  r,Ar  was  determined  by  the  sodar  dead  lone 

650  m.  For  each  strobe  pulse  starting  from  i«  the  first  iteration  cT  was  calculated  wi  ou 
turbulent  attenuation.  Then  the  first  iteration  of  the  outer  scale  of  turbulence  4°  was  calculated 
hv  Fn  (7  22)  Its  value  was  used  to  calculate  the  first  iterations  of  the  attenuation  coefficients  by 
Eos  (7  11  and  (7  12)  It  was  substituted  into  Eq.  (7.10)  to  calculate  the  first  iteration  of  the 
turbulent  attenuation  factor  which  was  further  substituted  into  Eq.  (7.1)  to  calculate  the  second 
iteration  cf2)(4>.  and  so  on.  The  iteration  process  terminated  when  the  number  of  iterations 
exceeded  100  or  when  the  next  iterations  started  to  diverge. 
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,vr  ?r  rhe  facsimilo  record  of  the  Zvuk~2  sodar  operated  on  September  29,  1996  from  14:55  to  16:47 
Tomsk  local  tile.  The  arrows  indicate  times  of  starting  of  three  subsequent  run;  then  data  were  processed 
with  the  use  of  the  proposed  iterative  algorithm. 

Figure  26  shows  the  facsimile  record  of  the  Zvuk-2  soda,  opcratecl  on  Septemto  29  19%  under 
conditions  of  developed  convection.  Intensive  thermal  plumes  up  to  120  200  m  can  be  seen,  in 
rowsTndicatc  times  of  starting  of  three  subsequent  lO-mln  runs;  the  eorrespondtng  vertical 

profiles  of  c]  are  shown  in  Fig.  27. 

For  comparison,  Fig.  28  shows  the  vertical  profiles  of  C\  reconstructed  from  the  data  of  acoustic 
sounding  in  the  morning  under  conditions  of  developing  convection.  The  structure  characteristics 
were  normalized  to  Cj  (*<>>  measured  at  the  minimum  sounding  altitude.  Their  numerical  values 

are  indicated  near  the  curves. 
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FIG.  27.  Normalized  vertical  profiles  of  c\  for 
three  subsequent  runs  indicated  in  Fig.  26. 


FIG.  28.  Normalized  vertical  profiles  of  C \  for 
three  subsequent  10 -min  runs  under  conditions  of 
de  ve l oping  con  vectio  n. 
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FIG.  30.  Vertical  profiles  of  the  outer  scale  of 
atmospheric  turbulence  for  three  subsequent.  10  min  runs 
on  September  29,  19%.  The  measurements  started  at 
10:00,  Tomsk  local  time,  under  conditions  of  developing 
convex  Hon.  The  corresponding  profiles  of  c\  are  shown 
in  Fig.  28. 
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The  solid  curves  are  for  the  first  iteration  and  the  dashed  curves  are  for  the  last  iteration  for  5% 
level  of  convergence  of  the  algorithm.  The  long-dashed  curve  illustrates  the  dependence 
C2T  (z)  ~  z~i/3  predicted  theoretically  for  convective  conditions.  All  curves  reach  their  minima  at 
2min.  Below  zmjn,  the  behavior  of  Cj  is  described  well  by  the  theoretically  predicted  dependence. 
Above  2min  c\  increases  with  the  altitude,  which  was  also  indicated  by  Neff  W.D.,  1975.  The 
asterisks  indicate  the  values  of  c\  07  m)  calculated  from  the  data  of  simultaneous  measurements 
with  micropulsation  sensors. 

Considering  £"4/3  dependence,  they  agree  well  with  the  sodar  data. 

Figures  29  and  30  show  the  vertical  profiles  of  L0  for  these  runs. 

The  solid  curves  show  the  first  iterations  and  the  dashed  curves  show  the  last  ones.  The  obtained 
profiles  were  approximated  by  the  nonlinear  least-squares  method  using  the  relation 


Lq(z)  «A  +  B  zc. 


(7.23) 


The  values  of  the  corresponding  coefficients  and  the  raw  fitting  errors  8  are  also  shown  in  the 
figures.  The  smooth  solid  curves  illustrate  this  empirical  dependence. 

As  can  be  seen  from  Figs.  29  and  30,  the  outer  scale  of  turbulence  increases  with  the  altitude  This 
agrees  with  the  existing  models  of  the  outer  scale  of  atmospheric  turbulence  reviewed  in  (Lukin 
V.P.  et  at.,  1997).  Its  near-ground  values  A  change  from  1.98  to  14.68.  According  to  Lukin  V.P. 
et  al  1997  the  effective  outer  scale  of  atmospheric  turbulence  determined  by  minimization  of  the 
integral  quadratic  discrepancy  of  the  structure  functions  of  the  optical  wave  phase  fluctuations 
corresponding  to  the  model  profile  of  the  outer  scale  L0(z)  and  the  constant  value  of  the  outer  scale 
r-  ik.s  between  0  68-34.7  m  for  the  lowest  intensity  of  the  atmospheric  turbulence  corresponding 
to  the  best  conditions  of  vision  through  the  atmosphere  and  between  1.31-55  4  m  for  the  highest 
intensity  of  the  atmospheric  turbulence  corresponding  to  the  worst  conditions  of  vision  through  the 
atmosphere.  It  should  be  noted  that  the  data  of  sodar  measurements  are  within  these  limits. 

Within  the  80-140  m  altitude  range  the  vertical  behavior  of  the  outer  scale  shown  in  Fig.  29b 
follows  the  dependence  L0(z)  =  2-fx  .  In  Figs.  29a  and  c,  L0  remains  practically  unchanged  at 
altitudes  up  to  150-200  m.  From  Fig.  30  it  is  seen  that  the  outer  scale  of  atmospheric  turbulence 
changes  quasiperiodically  with  a  period  of  10  min,  because  the  vertical  profiles  of  Lq  shown  in 
Figs  30a  and  c  are  very  similar.  That  is,  in  analogy  with  Lukin  V.P.,  1992,  we  suggest  that  10- 
min  intermittence  of  atmospheric  turbulence  (dramatic  changes  of  the  low-frequency  section  of  the 
turbulence  spectrum)  occurs  in  the  atmosphere. 

The  fact  that  in  the  first  run  under  conditions  of  developed  convection  (see  Fig.  29)  we  recorded 
the  most  strong  echo  signals  from  altitudes  up  to  550  m  has  engaged  our  attention;  at  that  time 
Ln  =  14.64  m  in  the  lower  part  of  the  atmospheric  boundary  layer  (ABL).  In  the  second  run,  the 
echo  signal  was  weak  and  Lq  =  8.74  m  in  the  lower  part  of  the  ABL  (see  Fig.  26).  In  the  thir 
run,  the  echo  signal  power  increased  and  Lq  =  9.57  m  m  the  lower  part  of  the  ABL. 


In  conclusion  it  also  should  be  noted  that  in  the  morning  under  conditions  of  developing  convection 
we  recorded  smaller  values  of  L0t  which  is  in  agreement  with  the  experimental  data  reported  in 
(Emaleev  O.N.  et  aL,  1979). 
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CONCLUSION 


The  effect  of  the  atmospheric  channel  on  the  near-ground  propagation  of  acoustic  waves  in  Je  atmosphere  has 
been  considered.  Numerically  shown  is  the  influence  of  the  meteorological  parameters  on  audibility  of  sound  ™ 
the  atmosphere.  The  structure  of  the  amplitude  and  phase  fluctuations  of  acoustic  waves  propagating  g 
turbulent  atmosphere  near  the  ground  has  been  analyzed.  Our  measurements  have  shown  a  significant 
dependence  of  the  temporal  characteristics  of  the  amplitude  fluctuations  on  the  angle  between  the  propaga 
path  and  the  direction  of  wind  velocity  vector.  Temporal  spectra  of  phase  difference  fluctuations  depended  on 
die  wind  velocity  much  weaker  than  in  case  of  line-of-sight  propagation.  The  structure  of  the  temP®r^ur®  “jj 
wind  velocity  fields  and  vertical  profiles  of  the  outer  scale  of  atmospheric  turbulence  have  been  experimentally 
investigated!  The  feasibility  for  model  descriptions  of  the  spectral  tensor  of  the  wind  velocity  field  fluctuations 
considered  as  homogeneous  anisotropic  turbulence  were  studied  based  on  the  experimental  data  obtained  by  us. 
For  the  first  time  the  vertical  profiles  of  the  outer  scale  of  atmospheric  turbulence  have  been  obtained  for  the 
atmospheric  boundary  layer.  Sodar-derived  outer  scale  values  agree  well  with  (he  estimates  of  the  effective  outer 
scale  of  turbulence  by  optical  methods. 
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